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Q . Abstract 

, We consider an immiscible two-phase flow in a heterogeneous one-dimensional porous medium. 

We suppose particularly that the capillary pressure field is discontinuous with respect to the space 
variable. The dependence of the capillary pressure with respect to the oil saturation is supposed 
to be weak, at least for saturations which are not too close to or 1. We study the asymptotic 
' behavior when the capillary pressure tends to a function which does not depend on the saturation. 

I In this paper, we show that if the capillary forces at the spacial discontinuities are oriented in the 

• , same direction that the gravity forces, or if the two phases move in the same direction, then the 

' saturation profile with capillary diffusion converges toward the unique optimal entropy solution to 

, the hyperbolic scalar conservation law with discontinuous flux functions. 
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1 Presentation of the problem 

The resolution of multi-phase flows in porous media are widely used in oil engineering to predict the 
^ motions of oil in subsoil. Their mathematical study is however difficult, then some physical assumptions 

■ have to be done, in order to get simpler problems (see e.g. [6, 11, 30]). A classical simplified model, so- 

^\ I called dead-oil approximation, consists in assuming that there is no gas, i.e. that the fluid is composed 

. of two immiscible and incompressible phases and in neglecting all the different chemical species. The 

oil-phase and the water-phase are then both made of only one component. 

- 1 — I 

^ . 1.1 the dead-oil problem in the one dimensional case 



Suppose that M represents a one dimensional homogeneous porous medium, with porosity (j) (which is 
supposed to be constant for the sake of simplicity). If u denotes the saturation of the water phase, and 
so (1 — u) the saturation of the oil phase thanks to the dead-oil approximation, writing the volume 
conservation of each phase leads to: 

ct)dtu + d,V^ = 0, (1) 
-(j>dtu + d^Vo = 0, (2) 

where Vo (resp. Vw) is the filtration speed of the oil phase (resp. water phase). Using the empirical 
diphasic Darcy law, we claim that 

Vf, = (O^p^ _ p^g) , p = o,w, (3) 

M/3 
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where K is the global permeability, only depending on the porous media, ^p,Pf:),pp are respectively 
the dynamical viscosity, the pressure and the density of the phase p, g represents the effect of gravity, 
fcr,/3 denotes the relative permeability of the phase /3. This last term comes from the interference of 
the two phases in the porous media. 

There exists Si, £ [0, 1) such that the function fc^,™ is non-decreasing, with kr_w{u) = if < w < 
< 1, and fcr.ti) is increasing on [s^^,!]. The function fc^^o is supposed to be non-increasing, with 

fcr,o(l) = 0. We suppose that there exists s* G (s*, 1] such that fcr,o(s) = for s G [s*, 1), and kr,o is 

decreasing on [0, s*). 

The pressures are supposed to be linked by the relation 

Pcap{u) ^P^-Po, (4) 

where Pcap is a smooth non-decreasing function called capillary pressure. 
Adding (1) and (2), and using (3) and (4) yields 

— {dxPi3 - Pl3g) = 0, 

and thus there exists 5, called total flow-rate, only depending on time, such that 

- E K'^idxP,~ppg)^q. (5) 

Using (3), (4) and (5), (1) can be rewritten 





Kdx ( , ) , ( ^ {dxPcapiu) - (p^ - Po).g) = 0. (6) 

Supposing that the total flow rate q does not depend on times, and after a convenient rescaling, 
equation (6) becomes 

ftw + 9x(/(M)-A(u)9,^(u)) =0, (7) 

where / is a Lipschitz continuous function, fulfilling /(O) = 0, /(I) = g, A is a nonnegative Lipschitz 
continuous functions, with A(0) = A(l) = 0, and tt is a non-decreasing function, also called capillary 
pressure. The effects of capillarity are often neglected, particularly in the case of reservoir simulation, 
and so (7) turns to a nonlinear hyperbolic equation called Bucklcy-Leverett equation, and we have to 
consider the initial-value problem 

dtu + dxf{u) = 0, 
u(0) = uo. 



1.2 discontinuous flux functions and optimal entropy solution 

We now consider heterogeneous one dimensional porous media, i.e. an apposition of several homo- 
geneous porous media with different physical properties. This leads to discontinuous functions with 
respect to the spatial variable. For the sake of simplicity, wc assume that the heterogeneous porous 
medium is made of only two homogeneous porous media represented by the open subsets Sli = R* and 
= Keeping the notations of (6), (j), K, kr,f3{u, •) and tt{u, •) are now discontinuous functions, i.e. 
piecewise constant functions, denoted 4>i, Ki,kr,p^i and tt^ in fi^. Thus the problem becomes 

dtu + dx{fi{u) - \i{u)dxTTi{u)) = 0, 

u(0) = wo, (8) 
+ transmission condition at a; = 0, 
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where fi are Lipschitz continuous functions on [0, 1], and can be decomposed in the following way: 

fi{u) ^ qri{u) + \i{u){py, - po)g, (9) 

where is a non-decreasing Lipschitz continuous function fulfilling ri[Q) ~ 0, ^■^(l) = 1, and A.^ is a 
non-negative Lipschitz continuous function fulfilling Ai(0) = 0, = 0. We stress here the fact that 

q and {pu, — Po)g neither depend on the subdomain i nor on time. 

We now have to give more details on this transmission conditions at x = 0. First neglect the 
effects of capillarity, so that (8) becomes the apposition of two Buckley-Leverett equations, linked by 
a transmission condition. 

dtu -t- dxfiiu) = 0, 

m(0) = uo, (10) 
+ transmission condition at a; = 0, 

We ask the conservation of mass at the interface between the two porous media, then we have to 
connect the flux. Denoting Ui the trace (if it exists) of u\a. on {x = 0}, this means that the following 
Rankine-Hugoniot condition has to be fulfilled: 

fliu,) ^ f2iu2). (11) 

Some assumptions has to be done on the flux functions fi in order to carry out the study. Firstly, 
we suppose that the total flow-rate g is a non-negative constant. Dealing with non-positive q is also 
possible, since it suffices to change x by —x and u by (1 — u). Secondly, we suppose that each /; has 
a simple dynamic on [0, 1]. More precisely, 

3bi E [0, 1) s.t. fi is decreasing on(0, &i) and increasing on {bi, 1). (12) 

With Assumption (12), we particularly ensure that 

9 = = max (/,;(s)). 

se[os] 

The physical meaning of (12) is that buoyancy works on the oil-phase in the sense of decaying x. 
The case hi = can also correspond to situations where the total flow rate q is sufficiently strong for 
ensuring that both phases always move in the same direction. Indeed, The oil-flux, given in fli by 
fi{u) has the same sign as the water- flux, given by g— fi(u). Note that the assumption on the dynamic 
on fi is often fulflUed by the physical models, as it is stressed in [1] (see also [29]). 

Thirdly, we assume 

/i and /2 are not linear on any non-degenerate interval of (0, 1). (13) 
This latter assumption allows us to claim, thanks to [38] (see also [44]) that a solution of 

dtu + d^fi{u) = 

satisfying the entropy inequalities in fli x (0, T): Vk £ [0, 1] 

dt\u-K\+d,{sign{u-K)if,{u)-f,iK)))<0 mV'{n,x[0,T)) (14) 
admits a strong trace Ui on {x = 0} x (0,T). 

Remark. 1.1 As it is proven in Section 3.1, if uq e BV{R), the solution u we consider has strong 
traces on the interface without assuming (13). The assumption (13) provides the existence of strong 
traces for initial data in L°° (R) . 

The problem (lO)-(ll) has been widely studied recently (see e.g. [1, 3, 5, 7, 8, 9, 22, 23, 24, 25, 31, 
32, 33, 34, 35, 36, 39, 41, 42, 43]). It has been particularly shown by Adimurthi, Mishra and Veerappa 
Gowda [2] that there are inflnitely many solutions satisfying (14). Additional entropy conditions has to 
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Figure 1: example of fi fulfilling (12) and (13) 



be considered at the interface {x ~ 0}. We refer to [2] and [14] for a detailed discussion on the possible 
choices of entropy conditions at the interface. According to Kaasschieter [33] and Adimurthi, Jaffrc and 
Veerappa Gowda [1] , the relevant entropy condition at the interface for two-phase flows with continuous 
capillary pressure field is the so-called optimal entropy condition introduced in [2] . Assuming that both 
fi are convex, the optimal entropy solution is characterized as follows: the discontinuity at the interface 
between ui and U2 can not be undercompressive: 

min{0,/((iii)}max{0,/^(?X2)} = 0. (15) 

Following the idea of Audusse and Perthanie [5] , the entropy condition at the interface can be derived 
by comparing the solution to steady states. Denoting by k{x) ~ Ki if a; e f^i, where /i(ki) = f 2(1^2)7 
then k has to satisfy the relation 

min{0,/KA£i)}max{0,/^(Ai2)} -0. 

We denote by Kopt{x) be the steady state corresponding to the optimal entropy connection appearing 
in the work of Adimurthi, Mishra and Veerappa Gowda [2], described on Figure 2 and defined as 
follows: 



if/i(foi) </2(&2),then 



if/i(?^i)>/2(62),then 



.optvx;-^ 6i=max{z.|/2(t^) = /i(6i)} if x > 0. > 

It is easy to check that in both cases, the function {x,t) '—^ /topt(a;) is a steady entropy solution in the 
sense of Definition 1.2. According to [2] and [14], the condition (15) can be replaced by the formulation: 

V^' e V+{R X [0,T)), 



\u ~ Koptldtip dxdt + / |mo - Kopt|V'(-,0) 
Jr 

y2 sigTi(u - Kopt) ifiiu) - /i(Kopt)) dxi'dxdt > 0. (18) 
Jo i^^Jn, 



This interface entropy condition does not require anymore that both fi are convex to provide a L^- 
contraction semi-group, as it will be stated in Theorem 1.3 and shown in [14]. 
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case (a) 



case (b) 



Figure 2: We exhibit particular steady states, which are so called optimal connections in the work of 
Adimurthi, Mishra and Veerappa Gowda [2]. The case (a) corresponds to the optimal connection (16) 
while the case (b) corresponds to the optimal connection (17). 



Definition 1.2 (entropy solution to (lO)-(ll)) Let uq € L°°(R) with 0<uq <1, and let T > 0. 

A function u is said to be an entropy solution to (lO)-(ll) if 

1. u€ L°°(]R X (0, T)), 0<u<l a.e.; 

2. for alii: eV+{Rx [0,T)), 



nudtipdxdt+ I uoip{-,0)dx+ ^ [ f 



f,iu)d,^dxdt = 0- (19) 



3. for i = 1,2, for all %p e V+ {Vl, x [0,T)), for all k G [0, 1], 



\u — Kldftp dxdt + / \uq — K\'ijj{-,0)dx 

T 

sigiliu- K){f,{u)- f^{K))d^^P dxdt >0; (20) 

J Hi 



4- the inequalities (18) hold. 



In the following theorem, we claim the existence and the uniqueness of the entropy solution to the 
problem (lO)-(ll). 

Theorem 1.3 (Existence and uniqueness of the entropy solution) Let uq G L°°{K.) with < 
Uq < 1, and let T > 0, then there exists a unique entropy solution to (lO)-(ll) in the sense of 
Definition L2. 

Furthermore, the function u can be supposed to belong to C{[0,T]]L\^J^)), and ifu,v are two entropy 
solutions associated to initial data Uq, Vq, then, for all R > 0, the following comparison principle holds: 

yte[o,T], 

rR rR+Ct 

/ {u{x,t) — v{x,t)y^dx < I {uo{x) — vo{x)y^dx, 

J-R. J-R-Ct 

where C = ma.Xi{Lip{fi)), with Lip{fi) sup^gjo |/j'(s)| ■ 

The uniqueness and the L^— stability of the entropy solution stated above can be seen as a straight- 
forward generalization of Theorem 3.1 in [14] to the case where g 7^ 0. The existence of such an 
entropy solution is provided in [1] by showing the convergence of the discrete solution corresponding 
to a Godunov-type scheme, while a modified Engquist-Osher scheme is considered in [14]. 
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In the particular case where fi{u) = kiu{l — u) for fci, k2 G M+, then it is shown in [41] that entropy 
solutions can be obtained as the limit for /i ^ and (5 ^ of the solutions u^-^ to the problem 

where is a smooth approximation of the piecewise constant function defined by k{x) = ki if .t e £7^. 

It has been proven in [33] that the entropy solution can also be obtained as limit for s tends to 
of weak solutions to regularized problems 

dtu' + 9, iMu') - eX,{u')d,7r,{u')) = 

under the assumption that 7ri(0) = 7r2(0) and 7ri(l) = 7i'2(l)- This latter assumptions is relaxed in this 
paper and in [16]. 

The Godunov-type scheme proposed by Adimurthi. Jaffrc and Veerappa Gowda [1] uses the fact 
that the fluxes are given by simple algebraic relations. Indeed, The flux G'int(wi, U2) = /i(wi) = f2{u2) 
at the interface corresponding to traces ui, U2 is given by: 

Gi„t(ui,U2) = max{Gi(ui,l), 62(0,^2)}, (21) 
where Gi is the Godunov solver corresponding to /,;, that is 

{min fi{s) if u < v, 
max ji[s) it u > V. 
sS lv,u] 

If the flux at the interface is given by 

Gint(^^l,W2) =/i(^/i) =Gi(wi,l), (22) 
then the restriction u\^^ of u to Qi is the unique entropy solution to 

dtu + d,fi{u) = inr!ix(0,r), 

u(0,-)=7 in(0,T), (23) 

u{-, 0) = uq in fli 

corresponding to 7 = 1. Recall that the trace on {x = 0} has to be understood in a weak sense (see 
[10, 37]). Since the solution to (23) depends in a non-decreasing way of the prescribed trace 7, we can 
claim that 

'^bi = sup^GL~((o,T);[o,i]) solution to (23)} . (24) 
Similarly, in the case where the flux at the interface is given by 

Gintiui,U2) - f2iu2) - 63(0,^3), (25) 

the restriction it|jj^ of u to ^2 is characterized by 

= inf76i-((o,r);[o,i]) l^' solution to (27)} , (26) 

where 

dtu + d,f2{u) = inr!2x(0,T), 

u(0,-) = 7 in(0,T), (27) 

ii(-, 0) = Wo in Q2 

Since the problem is conservative, in both cases the solution u is entropic in both fli x (0,T), i.e. it 
satisfies (14), and minimizes the flux through the interface {x = 0}. It is shown in [16] (see also [18]) 
that this characterization still holds, but that the different physical assumptions lead to the selection 
of a solution to (lO)-(ll) which is not the entropy solution in the sense of Definition 1.2. 
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1.3 heterogeneities involving discontinuous capillarities 



Let us now come back to the problem (8). Suppose for the sake of simpUcity that the functions tt^ 
are smooth and increasing on [0, 1], and that Xi{u) >OifO<u< 1. The problem is then a spatial 
coupling of two parabolic problems, and we will need to ask two transmission conditions: one for the 
trace, and one for the flux. Concerning the latter, the conservation of mass yields a relation analogous 
to (11), which can be written with rough notations: 

/i(mi) - Xi{ui)dxTriiui) = 72(^2) - X2iu2)dxTT2{u2). (28) 

Let us now focus on the trace condition at the interface. In the case of heterogeneous media, the 
capillary pressure can be discontinuous at the interface. Numerical schemes for simulating such flows 
has been proposed in [17, 26, 27, 28]. It has been shown independently in [15] and [20] (but see also 
[12] and [26]) that the connection of the capillary pressures Tri{ui) has to be done in a graphical sense, 
so that phenomena like oil trapping can appear. Thus we have to define the monotonous graphs tt^. 

{'K,{u) if < u < 1, 

{~oo,7Ti{Q)] if u = 0, 
[7ri(l),+oo) ifu=l. 

It is shown in [15] and [20] that a natural way to connect the capillary pressures on the interface 
consists in asking: 

^i(ui)n ^2(^2) 7^0. (29) 

In order to state a convenient definition for the solution of (8)-(28)-(29), we introduce the Kirchhoff 
transformation iy9i(u) = \i{s)Tr^{s)ds. 

Definition 1.4 (bounded flux solution) Let uq G L°°{M.), < < 1, and let T > 0. A function 
u is said to be a bounded flux solution to (8)-(28)-(29) if it fulfills: 

1. ue L°°(]R X (0, T)), < It < 1 a.e., 

2. d^ip,{u) e L°°(r!, X (o,r)), 

S. ^ri(ui)n7ir2(u2) 7^0 /or a.e. te{Q,T), 

4. yi' G v{R X [o,r)), 

/ / u{x,t)dtip{x,t)dxdt + / uo{x)ifi{x,0)dx 
Jo Jr Jr 

+ [ J2 [ {h{u)[x,t)~d^v^{u)[x,t))d^^{x,t)dxdt = Q. (30) 



The bounded fiux solution are so called since the point 2 of Definition 1.4 ensures that the flux 
fi{u) — dxfi{u) remains uniformly bounded. Such a condition will require assumptions on the initial 
data Mo, as it will be stated in the following theorem. 

Theorem 1.5 (existence of a bounded flux solution) Let /i,/2 be Lipschitz continuous func- 
tions, and (pi,ip2 be increasing Lipschitz continuous functions. Let uq G L^{M.), < uq < 1 fulfilling 
dxipi{uo) G L°°{ni), and 7ri(Mo^i) n7r2(Mo,2) 7^ 0, where uo,i denotes the trace on {x = 0} o/uqiq. . 
Then there exists a bounded flux solution. Moreover, u belongs to C([0, T]; L^(]R)). 

The first part of this theorem is a straightforward adaptation to the case of unbounded domains and 
non- monotonous fi of a result from [17] and [20] (see also [40]). This is based on a maximum principle 
on the fluxes (/i(u) — dx^Pi{u)). This particularly yields: 

\\fi{u) - 9r'^t(w)|lL~(n,x(o,T)) < max (||/j(mo) - 92;<^j(uo)I!l~(Oj)) • (31) 

3 — li-^ 
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If 1*0 e L'^i^), then choosing 7/' = mm(l, (1, i?-|a;|)+) and letting i? tend to 00 gives it € i°°((0, T); ^^(R)). 
Moreover, thanks to [19], u can be supposed to belong to C{[Q,T]; Ll^^{M.)). Then u belongs to 

C([0,T];Li(K))- 

The choice of bounded flux solutions instead of more classical weak solution with dx^i{u) only 
belonging to i^((0, T); Lf^^{ni)) has been motivated by the fact that it provides a comparison principle. 

Proposition 1.6 Let u,v be two bounded flux solutions in the sense of Definition I.4 associated to 
initial data uq^vq- Then, for all ip G P^(M x [0, T)), 



{u{x,t) — v{x,t)) dti'{x,t)dxdt + / {uo{x) — vo{x)) ip{x,0)dx 
I Js. 

+ 1 X! / sign^{u{x,t)~v{x,t)){fi{u){x,t)- f,{v){x,t))dx'ilj{x,t)dxdt 

/ V/ d,{^,iu){x,t)-^,{v){x,t))^d,ij{x,t)dxdt>0. (32) 



i=l,2 



This proposition is not sufficient to claim the uniqueness, but it will be very useful in the sequel. In 
order to obtain a uniqueness result, we have to ask furthermore that the initial data belongs to L^(M). 

Theorem 1.7 (uniqueness of bounded flux solution) Let uq G L^{M.), < mq < 1 a.e., with 
dxy^iiuo) € L°°{fli) and 7ri(ito,i) H 7r2(uo,2) 7^ 0. Then there exists a unique bounded flux solution 
u E C{[0,T]; L^{W)) in the sense of definition L4-. 



This theorem is a straightforward consequence Proposition 1.6. Indeed, choosing -0 = min(l, (1, R — 
|a;|)+) in (32), and letting R tend to +00 gives the comparison principle: g [0,T], 

/ {u{x,t) — v{x,t))^ dx < / (uoix) — vo{x))^ dx. (33) 
Jr Jr 

The uniqueness result follows. 



1.4 capillary pressure independent of the saturation 

In some cases, the dependence of the capillary pressure tt^ with respect to the saturation seems to 
be weak, and some numerical simulation consider capillary pressures only depending on the porous 
medium, but not on the saturation. More precisely, we aim to consider graphs of capillary pressure on 
the form 

( P, if < u < 1, 

^^(^t) = <^ i-oo,P,] if u = 0, (34) 
[ [P, ,+oo) if u = 1, 

so that the capillary pressure would roughly speaking not depend on u. 

If one considers an interface {x = 0} between two H^, where the TTi are on the form (34), we can 
give an orientation to the interface: the interface is said to be positively oriented if Pi > P2, and 
negatively oriented if Pi < P2- A positively oriented interface involve positive capillary forces, and a 
negatively oriented involve positive capillary forces. The gravity effects are also oriented by the sign 
of {pw ~ Po)g in (9). We have to make the assumption that 

"either the gravity effects and the interface are oriented in the same way, or the convective effects are 

larger than the gravity effects. " 

Since we have supposed that gravity works in the sense of decaying x, we assume in the sequel that 

Pi < P2. (35) 
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We build a family of approximate problems {T"^) taking into account the capillary pressure: one 
suppose that 7r|(ti) = Pi + eu, where Pi is a constant depending only on the homogeneous subdomain 
fli. In fact, any nf converging uniformly to Pt on [0, 1] and such that u i-^ JJ' Ai(s) (nf)' {s)ds converges 
uniformly toward would fit. 

Up to a smoothing of the initial data, we obtain a resulting sequence (u'^)^ of bounded flux solutions 
for a problem of type (8)-(28)-(29). We will show that under Assumptions (12)-(13), this sequence 
tends almost everywhere to the unique entropy solution to (lO)-(ll). 

This result has to be compared to the one presented in the associated paper [16], where it is shown 
that if Assumption (35) does not hold, non-classical shock can occur at the interface, representing 
oil-trapping. 

1.5 organization of the paper 

The paper is organized as follow: section 2 is devoted to the study of the approximate problem (V^). 
We first smooth the initial data in a convenient way, and then we give a L^((0, T); i/^(rii))-estimate 
on the approximate solutions that shows in particular that if the approximate solution u*^ converges 
almost every where towards a function u, then u satisfies the points 1,2 and 3 of Definition 1.2. In 
order to prove that converges almost everywhere, we derive a family iJT^-estimates. In order to 
check that the last point of Definition 1.2 is fulfilled by the limit u of the approximate solutions (u^)^ 

2 The approximate problems 

In this section we will define the approximate problem {T"^), and its solution u^. We will state a 
L^((0, T); iJj^Q^(ni))-estimate and a family of BF-estimates, which will be the key points of the proof 
of convergence of u'^ toward a weak solution of the problem {T"^). 

In order to recover a family of entropy inequalities, we will build some steady solutions to the 
problem (V^), and study their limit as e — *■ 0. This last point will require strongly Assumption (12). 

2.1 smoothing the initial data 

As it has already been stressed in Theorem 1.5, we need to assume some regularity on the initial data 
to ensure the existence of a bounded flux solution to problems of the type (8)-(28)-(29). 
Let Mo belong to L°°(M), with < wo < 1, we will build a family {uq)^ of convenient approximate 
initial data. 

Lemma 2.1 Let uq G i°°(M), < 5: 1; then there exists a family {uq)^ of approximate initial data 
such that: 

• w5 e C,°°(R*), < u5 < 1, 

• Uq —^ Uq a.e. in R, ||e9a;itQ||oo — > as e — > 0, and Wsd^UQWac < 1 for all e > 0, 

• IfuQ e BV{R), then !|a,ug|Ui(R) < TV{un)+i. 

Proof: Let a > 0, and let pa be a mollifier with support in {—a, a). Let = (uq o Xa<\x\<i/a) * 
Pa, then it is clear that G C^(M*), and that — > uq a.e. in M as a 0. Choosing e = 

min ^a, , a-nd Uq ~ v" ends the proof of Lemma 2.1. □ 

2.2 the problem (V^) 

Let Pi, P2 G M, we define the functions nf by 7r|(ti) ^ Pi + eu, and 

{Pi+eu if < u < 1, 

(—00, Pi] if u = 0, 
[P,+£,+oo) ifu=l. 
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If e is small, the intersection of the ranges of the functions 7rf and tt^ is empty, and then the graphical 




Figure 3: capillary pressures graphs nf. 



relation 7ri(wi) n tt2{u2) 7^ connecting the capillary pressures at the interface becomes 

(1 - Ui)u2 = 0. 

Let < uo < 1; and let (uq)^ be built as in Lemma 2.1. let (pi{u) ~ Xi{s)ds. The approximate 
problem is: find s.t. 

dtu' + iMu') - ed^ip,{u')) =0 in r!, X (0,T), 

^JK)(o-,t)n^?^(^z^)(o+,<)^0 in(o,r), 

h{u'){0- ,t) ~ ed,ipi{u')(0- ,t) ^ f2{ul){Q+ ,t) ~ ed,^2{u'){Q+ ,t) in(0,r), ^' > 

w^(0) = u% in M. 



This problems is of type (8)- (28)- (29), then the notion of bounded flux solution is a good frame to 
solve it. 

Definition 2.2 (solution to (V^y) A junction is said to be a (bounded flux) solution to (V^ ) if 
it fulfills 

1. e L°°(]R X (0,T)), < < 1 a.e., 

2. d^^,{u^)eL^{Sl,x{0,T)), 

3. VV' e 2?(M X [o,r)), 

T 

u'^(x,t)dt'il'ix,t)dxdt + I ul{x)ip{x,0)dx 



1=1,2 



^ / {f,{u^){x,t)-ed^^,{u'){x,t))d,il:{x,t)dxdt = 0. (36) 



We can use Theorem 1.5 to claim that there exists a family (m^)^ of bounded flux solution to {V^) in 
the sense of Definition 2.2. Moreover, this family of solution fulfills, thanks to (31) and Lemma 2.1: 
for all e > 0, 

\\edxipi{u'^)\\oo < (max(Lip(.^i)) + max||/i||ioo(o,i)). (37) 
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Since Uq belongs to L^(R), the solution u'^ is furthermore unique in C([0, T], L^(M)) thanks to Theo- 
rem 1.7. 

2.3 a L2((0,T);/f/„^(ni))-estimate 

All this subsection is devoted to prove the following estimate. 

Proposition 2.3 Let K be a compact subset of Qi, and let u'^ be a solution of (V^ ) in the sense of 
Definition 2.2, then there exists C depending only on fi,K,T (and not on e) such that 

Vs\\Vi{u'')\\L^-((0,T)-m(K)) < C. 
Particularly, this implies that edx^i{u^) — > a.e. in fi^ x (0, T) as e — > 0. 

Proof: We fix e > 0. Since the functions ip^^ are not Lipschitz continuous, the problem (V^) is not 
strictly parabolic, and the function u"^ is not a strong solution. In order to get more regularity on the 
approximate solution, we regularize the problem by adding an additional viscosity 1/n {n > 1), so 
that the so built approximate solution is regular enough to perform the calculation below. 

Let n > 1, and let (pi.n{u) = (pi{u) + u/n, and let be a bounded flux solution of (V^) with 
ipi^n instead of ipi. From (37), we know that dxfi,n(u'^n) is uniformly bounded in L°°{Qi x (0,T)), and 
since ip~^^ is a Lipschitz continuous function, one has dxU,-^ € L°°{fli x (0,T)). The following weak 
formulation holds: VV^ e V{n x [0,r)), 

T 

u'^{x,t)dtip{x,t)dxdt + I UQ{x)^jj{x,0)dx 



J2 / {f,{O{x,t)~edx^^,n{u'J{x,t))dxijix,t)dxdt = 0. (38) 



Let (a, 6) C r^i. Let ( e I?+((a, b)), we deduce from (38) that 



(at<|<C')=/ / /,«)a,«C')da;dt-e / / ^xip^,ni<)^A<(:^)dxdt, (39) 



where ^ • | • y is the duality bracket between L'^{{0,T); H ^{a,b)) and L^((0, T); iJg (a, 6)). Since ipi^n 
is a Lipschitz continuous function with (||Aij|oo + 1/"-) as Lipschitz constant, one has 



Tpb ^ pT i-b 



^xip^,n{O^AOC''dxdt>———— / idx^,.„{u'jf edxdt. (40) 



'OJ a 

Let $i be a primitive of fi, then: 



I I f.{u'Jdx{ut,e)dxdt =11 dxM<)C^dxdt+ f f f,{ul^)u'M''dxdt 

[/,«)< -$,«)]a,C'da;rfi. (41) 



Admit for the moment Lemma 2.4, stated and proven below. We deduce from (39), (40), (41) and 
Lemma 2.4 that 

T rb 

{dxtpi,niut,)f C'^dxdt 



< h\C\\l\b -a\+ [ [ |< (/.«) - sdx^^AO) - M<)\ \d.e\dxdt. 



2' 



J a 
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Using now the fact that u^^ is a bounded flux solution, we deduce from (31) that [un{fi{un) —edx'Pi,n{'Un))~ 
^j(^n)] uniformly bounded independently of e and n, and so there exists C only depending on fi, 
\b — a\, uq and such that: 

{d^^.AOf Cdxdt < C {\\d,e\\LHio,Ty,Mm) + IICIIL) ■ 

2 

This estimate still holds for ({x, t) = X{a,b)(x), for all (a, 6) G fi^ , so we obtain 

(a,v5,,„«))' dxdt < C{2T + 1). (42) 

Classical compactness arguments provide the convergence, up to a subsequence, of (w^)„ to a solution 
of {V^) in LP(M. x (0, T)), 1 < p < +oo. This ensures particularly that, up to a subsequence, 

lim < = u" a.e. in M x (0,T). 

n — >-\-oo 

Taking the limit w.r.t. n in (42) yields 

i-T rb 

{dxipiiu"))^ dxdt < C{2T+1). 



OJa 



□ 



Lemma 2.4 Let be an approximate solution of (P^ ) with ipi^n instead of ipi, and let G I?^((a, b)), 
then 

{dtut,\uie)>-^\\C\\l\b-a\. 

Proof: Since «C) e L^iiO^T); H^{a,b)), and St «C) S L^iiO,T); H-\a,b)), and so, up to a 
negligible set, «C) e C{[0,T]; L\a,b)), and 



= <(x, T)2c'(2:)dx - 1 ^.o(:r)2c2(a;)dx 

> -^IICIlLl^-al- 



□ 



2.4 the SF- estimates 

In this section we suppose that uq G BV{M.). In order to avoid heavy notations which would not lead 
to a good comprehension of the problem, the following proof will be formal. To establish the following 
estimates in a rigorous frame, one can introduce of a thin layer {—rj, rf) on which the pressure variates 
smoothly to replace the interface, and add some additional viscosity to obtain smooth strong solutions 
to the problem. This regularization of the problem has been performed in [12] and in [20]. 

For a, 6 G [0, 1], we denote by 

Fi(a, b) = sign{a - b){f^{a) - fi{b)). 



Lemma 2.5 There exists C depending only on fi, T , uq such that 

|9(i^i(u^, k)|^^(j^^^(p j,^) < C. 
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Proof: Suppose in the sequel that is a strong solution, i.e. 

holds point-wise in fli x (0, T). Let h > 0, and let t € {0,T — h). Comparing u'^{-,- + h) and u'^ with 
(33) yieds 

\u'^ (x , t + h) — (x , t)\dx < j \u^{x,h) — UQ{x)\dx. 



Dividing by h and letting h tend to 0, one can claim using the fact that is supposed to be a strong 
solution 

/ \dtu'{x,t)\dx ^ I \dA.h{u'){x,t) - ed^^,{u'){x,t)]\dx 



1=1,2 



< 



\dtul{x)\dx 



- X! / \dx[fi{ul){x) - £^^^p^{ul){x)]\dx 

1=1,2 



Lemma 2.1 then ensures that there exists C not depending on e such that 



\dtu%x,t)\dxdt = / V / |9:.[/,K)(x,t)-£a,¥'«("')(^,t)]Mxdi < C. (43) 



Thanks to the regularity of /i, this particularly ensures that, if we denote by M^iyti x (0,T)) the set 
of the bounded Radon measure on O,; x (0, T), i.e. the dual space of Cd^i x [0, T), M) with the uniform 
norm, we obtain: Vk € [0, 1], 



oo • 



Lemma 2.6 There exists C depending only on uq, fi and T such that 



□ 



[Fi{u^,ti) - edx\'Pi[u^) - <yJj(K)|) 



< C. 



Mt{niX{0,T)) 

Proof: It follows from the work of Carillo (see e.g. [21]) that for aU k e [0, 1], for aU V' e X'+(a x [0, T)), 

l-u*^ - nldtipdxdt + - n\ij{0)dx 

Jn, 

Fi{u'',K) - £dx\^Pi{u'') - Lpi{K)\\dxil^dxdt > 0. 



T 

J Hi 



T 

J Hi 



(44) 



Let 77 > 0, we denote by uJrj{x) = (1 — \x\/ri)'^ ^ and suppose now that ij} belongs to V'^ijli x [0,T)), i.e. 
ij} does not vanish on the interface {x = 0}. Estimate (44) still holds when we consider tpjj = ip{l — LOn) 
as test function. 



\u'^ ~ k\{1 — uj,-i)dt'il'dxdt + / \u'^ — k\'\Jj{0)[1 — ujr^)dx 
Jni J^i 



Jiii 

T 



Fi{u'' ,k) - edx\'^i{u'') - <(5j(k)| ) (1 - u)n)dxi'dxdt 



> 



Jn 



Fi{u^, k) - edxlipiiu") - (fi{Hi)\]ipdxUj,jdxdt. 



(45) 
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The fact that the flux induced by m*^ is uniformly bounded w.r.t. e thanks to (37) implies that there 
exists C depending only on uq, fi such that 



Then we obtain the following estimate on the right-hand-side in (45): 



(46) 



Jfl 



dxdt 



Letting i] tend to in inequality (45) gives: V?/' G V+i^i x [0,T)), Vk G [0, 1], 



\u'^ — K\dtipdxdt + — K\ip{0)dx 

JUi JUi 



T 

Jfli 



(i^,(u^K) - ed^\ip^iu') ~ ^,{K)\^d,^pdxdt > -CTU'\ 



(47) 



We introduce now a monotonous function Xip G ^^{Qi) equal to 1 on the support of il>{-,t) for all 
t e [0,T] (so that Wd^xALHUi) = 1), then ||i/'||ooXv ^ '^^ Choosing ||'(/'||ooXv " V' and ||V'||ooXV + as 
test function in (47) yields, using (46) once again 



9t|u^ - k\- dx{Fi{u^,K) - £dx\ipi{u'') - ^Piin)]^ 



< 2CT. 



(48) 



A^6(n.x(0,T)) 

Thanks to (43), there exists C depending only on ito, fi and T such that 

\dt\u^~^\\M,iiMMo,T))<C'- (49) 
Lemma 2.6 is so a consequence of (48) and (49). □ 

and let K = [a, 6] C fi^. We introduce z]^{x,t) defined on the 



Proposition 2.7 Let ua e 

whole space given by: 



F,iu^,K)ix,t) if{x,t)eKx{0,T), 
otherwise. 



There exists C depending only on uq^ fi^T^ I\ and o, unifonnly bounded function tj^^ i^, with rK,t,{£) 
tends uniformly to with respect to k as e —* 0, such that, for all (^, h) G R^, 



I^kA^ + ^,t + h)- zl,Jx, t)\ dxdt < C(|e| + \h\) + rK,«(e). 



Proof: Let h G M., one has 



\zK,Ki^^^ + ^) - ZKAx,t)\dxdt < |9tz^„|^^^jj2) 1^1 
< (|at^^,(^/^«)|^^(^^,(o_^))+2|lF,(^^^K)|U.^(T + fe-a)) \h\ 
It follows from Lemma 2.5 that there exists Ci depending only on uq, fi, K, T, ipi such that 

/ / \z%L K,i^' t + h) — zl^ ^(x, t)\ dxdt < Ci\h\. 

We also define 



(50) 



ed,]^,{u^){x,t)-^d.'i)\ ii{x,t)eKx (o,r), 
otherwise. 
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Proposition 2.3 ensures that qj^ ^{x,t) converges to almost everywhere in R^, and the estimate (46) 
ensures us that g|f ^(^i stays uniformly bounded in L°°(R^) with respect to e. 
Let ^ € K, then we have: 



< 



Using (46) and Lemma 2.6 in (51) yields that there exists C2 depending only on uq, fi, K,T, ipi such 
that 



\zKjx + ^,t) - (?^,«(cc + e,t) - {zl,Jx,t) - qjijx,t)) I dxdt <C2\^\. 
This particularly ensures that: 

\zKAx + ^,t)-z'^Jx,t)\dxdt<C2m+2\\qj,^,ix,t)\\m^2y (52) 

By choosing C = max(Ci, C2), one we deduce from (50) and (52) that 

\z'kJx + ^,t + h)- z'j,Jx, t)\ dxdt < Ci\h\ + lel) + 2\\ql,Jx, t)|Ui(K2). 



We conclude the proof of Proposition 2.7 by checking that \\qj^ ^(a;, t)||Li(i{2>) converges uniformly to 
with respect to k as e tends to 0. □ 



2.5 some steady solutions 

The spatial discontinuities of the saturation and capillary pressure allow us to consider Kruzkov en- 
tropies \u — k\ only for non-negative test functions vanishing on the interface, i.e. in 'D'^{R* x [0,T)). 
This is not enough to obtain the convergence of u"^ toward an entropy solution u in the sense of 
Definition 1.2, and a relation has also to be derived at of the interface. 

In order to deal with general test functions belonging to V'^ (M x [0, T)), we will so have to introduce 
some approximate Kruzkov entropies | • — k'^(x)|, where k*^ are steady solutions of (V^). Letting 
£ — > 0, those steady solutions converge to piecewise constant function defined below. The functions 
I ■—k^{x) \ correspond to the so called partially adapted entropies introduced by Audusse and Perthame 
[5] . We will then be able to compare the limit u of approximate solutions to this limit , and then 
to prove that u is the unique entropy solution. 

The building of convenient strongly uses Assumption (12). It will be shown in [16] that if (12) 
fails, non classical shocks can occur at the interface, and the limit u of the approximate solutions u'^ 
is thus not an entropy solution. 

Recall that q > 0, Pi < P2, and suppose that < e < P2 — Pi- Some simple adaptations can be 
done to cover the case q < 0. The transmission condition (29) can be summarized as follow: cither 
Til = 1, or U2 — 0. 

We have to introduce the following sets: 

£1 = {ki /3k2 with /i(ki) = /2(k2)}, 

£2 = {k2 with /i(ki) = /2(k2)}. 

It follows from Assumption (12) that either £1 ~ [0,1], or £2 ~ [0,1], and so we are ensured that 
K G [0,1] belongs either to £1, or to £2 (or of course to both). Check also that for all z e (0, g], it 



15 



follows from (12) that there exists a unique Ki(z) such that fi{hii{z)) ~ z. On the contrary, if z e £i 
with z < 0, 2; has two antecedents through Let k belong to one denotes 

n\ = max{i/ | U{v) = (53) 

and 



K?- = mm 



I h{y) = /,(^)}. (54) 

Definition 2.8 (reachable steady state) A function k{x) is said to be a reachable steady state if 
there exists a steady entropy solution k^{x) to the problem {V^) in the sense of Definition 2.2 converging 
to k{x) in L^(R) as £ tends to 0. 

This section is devoted to establish the following proposition, that exhibits all the reachable steady 
states. In particular, all the steady states that are not undcrcompressiblc are reachable. 

Proposition 2.9 For all k €z £j, there exists a family of steady solutions (k^)^ to the problem (V^ ) 
such that 

k'^ — > k-' a.e. in fti, (55) 

where (x) can be chosen between: 

i) (x) — K^ if x < and x > 0; 
ii) k^ (x) — K^ if X < and K2 if x > 0; 
Hi) k^ (x) = k{ if X < and Kg x > 0. 
In particular, Kopt defined in (16)-(17) is a reachable steady state. 

Proof: Let k Cz £j. 

• If K = 1, the three limits k^ are identically equal to 1, which is a steady solution fulfilling 
Proposition 2.9. 

• We suppose now that k < 1, and fjin) > 0. Even in this case, the three reachable limit are 
the same. Thus we only have to build one sequence of converging steady solutions. Let y be a 
solution of: 

-^^i{y) = fjii^) - hiy), forx>o, ^^^^ 

y(0) = 1. 

The solution y{x) converges to k{ as a; +00. The family (k"^)^ defined by: Ve 

K%x) = l ^2 0r K^2 ij^>0> (57) 
^ ' { y{-x/e) if a; < 0. ^ ^ 

fulfills so the conclusion of Proposition 2.9. 

• Suppose now fj^n) ~ 0. The solutions k'{x) built with (56)-(57) converges toward the two 
reachable steady states i) and ii). One can also choose k^{x) = 0, which is of course a steady 
solution. 

• It remains the case /j(k) < 0. The solutions K^(a;) built with (56)-(57) still converges toward 
the two reachable steady states i) and ii). 
Let m; be a solution of 

-^ip2{w) = f2iw) - fj{K), for .T > 0, 

u;(0) = 0. 
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and 

X _ J k{ if X < 0, 

w{x/e) ifx>0, 

then converges toward the third reachable steady state. 

□ 

3 Convergence toward the entropy solution 

3.1 almost everywhere convergence 

We state now a lemma which is an adaptation of Helly's selection theorem criterion, which states that 
if (v^);, is a family of measurable functions on an open subset U of R*' (fc > 1), uniformly bounded 
in L'^JU) and BViociU)i one can extract a subfamily still denoted by (w^)^ that converges almost 
everywhere in U, and the limit belongs to BViociU)- 

Lemma 3.1 Let U he an open subset of M.^ (k > 1). Let (w^)^ be a family of functions uniformly 
bounded in L°°{IA) with respect to e. Let K be a compact subset of 14. For C € M.^ , we denote by 

K(; = {xeK \x + CeK }. 

One assumes that there exists C > depending only on K ( and thus not on e) and a function r 
fulfilling linig^o = , such that for all C G R'"', 



\v^-{x + C) - v%x)\dx < C\C\ + r(e). (58) 



Then, there exists a sequence (£«)„ tending to 0, and v e BVioc{l^) such that 

lim u'" =v a.e. inlA. 

n — *+oo 

Proof: Let K he a. compact subset of U. Estimate (58) says, roughly speaking, that is almost a 
Sy-function on i.e. as close to BV{K) as wanted, provided that e is supposed to be small enough. 
So we will build a family [w^)^ of BV^-functions, which will be close to the family (w'^)^, at least for 
small e, and we will show that (w"^)^ admits an adherence value v in BVlI'C) for the L^(/C)-topology, 
and that this v is also an adherence value of {v'^)^- Another proof for the a.e. convergence toward 
a function v can be derived directly from Kolmogorov compactness criterion (see e.g. [13]). But the 
advantage of the following method is that it provides directly some regularity on the limit v G BViocipi). 

Let (p'^)j be a sequence of mollificrs, i.e. smooth, non negative and compactly supported functions 
with support included in the ball of center and radius e, and fulfilling |lp'^||Li(Rfc) = 1 for all e > 0. 
We define the smooth functions 

where v^{x) = v^{x) if x Cz U and v^{x) = if a: e W^. 
Thanks to the regularity of w^. 



\w 



-■{x + 0-w'{^)\dx<\\Vw^\^^^,^j,)).\C\. (59) 



Suppose that e < d{K,dU) (with the convention d{K,(d) = oo). Thanks to (58), we have also 



\w'{x + 0~w'{x)\dx < / \v'{x + C)-v'ix)\dx <C\C\+r{e), (60) 

where 

{K(; + e} = {xeU \ d{x, K) < e}. 
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Since r(e) tends to as e ^ 0, this particularly ensures 

limsupllVw^lUwRfc))!. < C. 

The family [w'^)^ is thus bounded in BV{Kc^) in the neighborhood of e = 0, and thus, thanks to Helly's 
selection criterion, there exist v G BV{Kq), and (£«)„ tending to such that 

w'^" — > V a.e. in K as n ~r oo. 

Furthermore, for all n G N*, 

llw''" - w^'iL^if) ^ \v^" {x - y) - v^" {x)\p''" {y)dydx 

< Ce„ + r(£„). 

This ensures that v^" tends also almost everywhere toward w as n tends to +oo. □ 
Lemma 3.1 will be used to prove the following convergence assertion. 

Proposition 3.2 Suppose that uq G BV{M.), and let he a solution to (V^ ). Up to an extraction, 
there exists u G L°°(M x (0, T)), < u < 1 a.e. such that 

u" a.e. m M X (0,r). 

Furthermore, there exists mi,M2 G L°°[Q,T), such that 

1 F f° 

lim - / / \u(x,t) — ui(t)\dxdt = 0, 
''-0 V Jo 

1 F r 

lim - / / \u{x,t) - U2{t)\dxdt = 0. 

Proof: Let K be a compact subset of 
We define the function H, : [0, 1] M by 

H,{u)^ [ {F,{u,a)- h{a))da, (61) 

so that, thanks to Proposition 2.7, there exists C depending on ug, T, K, and a function r tending 
to as £ tends to such that for all ^ G R, /i G (0, T), 

/ \H,{u'){x + ^,t + h)-H,{u'){x,t)\dxdt<C{\£,\ + \h\) + r{e). (62) 
Jk 

An integration by parts in (61) yields: Vu G [0, 1] 



H,{u) = - I {a-hi)d„{F,{u,a)- h{a))da + {2h-l)f,{u) 

{a-h)fi{a)da + {2h~\)U{u) (63) 



where, thanks to (12), fi is decreasing on [0,6^] and increasing on [5^, 1]. Using Proposition 2.7 with 
K = 0, 

/ |/.(^^^)(x + e,t + M-/.K)(s,<)M^d^<C(|?| + |/i|) + r(£), (64) 

JK 

where C and r have been updated. Denoting by 

Ai{u)^ I (cr - 6,)/;(a)dcr, 
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we obtain from (63) and (64) 

T-h !■ 

/ \A,{u'){x + tt + h)- A,{u^){x,t)\dxdt <C(\i\ + \h\)+r{e), (65) 
Jk 

with a new update for C and r. Thus wc deduce from Proposition 3.1 that, up to an extraction, 
Aj(it^) converges almost everywhere toward Ai e BV{K x (0,T)). It follows from (12) that Ai is an 
increasing function, and so we obtain the convergence almost everywhere in K x (0, T) of u'^ toward a 
measurable function v. 

Since for all e > 0, < < 1, there exists u G L°°(R x (0, T)), < u < 1 such that converges 
to u in the L°°(]R x (0, r))-*-weak sense, we thus have, up to an extraction, 

u" a.e. in K x (0, T). (66) 

Since (66) holds for any compact subset K of Qi, for i = 1, 2, we can claim that, up to an extraction, 

u" a.e. in M X (0,T). 

Moreover, since Ai{u) belongs to BV{M. x (0,T)), we can claim that Ai{u) admits a strong trace on 
{x = 0} X (0, r) (see for instance [4]). Using once again the fact that A~^ is a continuous function, 
we can claim that u admits also a strong trace on each side of the interface. □ 

3.2 convergence toward the entropy solution 

In this section, it is proven that the limit value u for the family (u^)^ exhibited previously is the 
entropy solution described in Definition 1.2. 

Proposition 3.3 Let Uq G BV{R), < Uq < 1 a.e., and let (ug)g a family of approximation of Uq 
given by Lemma 2.1, and let (u^)^ be the induced sequence of bounded flux solution of (8)-(28)-(29). 
Then under Assumption (12) 

lim u' in LPJR x [0, T]), Vp e [1, oo) 

e— >0 

where u is the unique entropy solution to (lO)-(ll) associated to initial data uq. 

Proof: Thanks to Proposition 3.2, we can suppose that there exists u G L°°(R x (0, T)) such that, up 
to a subsequence, 

u' u in L,oc(M X (0, T)) as £ ^ 0, (67) 

then 

f.{u') ^ f^{u) in Ll^in, X (0,T)) as £ 0. 
Furthermore, thanks to Proposition 2.3, one has 

ed,,^,{u') ^ in Ll^{n, x (0,T)) as £ 0. (68) 

As a consequence, letting e tend to in (36) provides (recall that Mq tends to uq in L\^^{M.)) that u is 
a weak solution to (lO)-(ll), i.e. that it satisfies (19). 

Since for £ > 0, {eLpi) ^ is a continuous hmction, it follows from the work of Carrillo [21] that it 
fulfills the following entropy inequalities: Vk G [0, 1], V?/' G 'D+{Vli x [0,T)) 

T , 

— hi\dtip dxdt + I \uq — {■,{)) dx 
/o JUi JQi 

+ f f {F,{u',n)-ed^\ip,{u')-ip,{K)\)d.^i:dxdt>0. (69) 

Jo J Hi 
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The assertion (67) particularly yields that for i = 1, 2 and for all k G [0, 1], 

Fi{u'',K) Fi{u,K) in Ll^{R x (0,T)) as e ^ 0. (70) 
On the other hand, it follows from (68) that 

ed,\^,{u') - ^,{k)\ ^ in Ll^{R x (0,r)) as e ^ 0. (71) 

Taking (67), (70) and (71) into account in (69) provides the inequality (20). 

The last point remaining to check is that the interface entropy condition (18) holds. Let (k'^)^ be 
a family of steady states to {T"^) converging in Lj^^{M.) as e tends to to Kopt defined in (16)-(17). 
Recall that such a family exists thanks to Proposition 2.9. Then for all fixed e > 0, is a steady 
bounded flux solution. Hence, it follows from (32) that 

T f 

lu" - K''\^t^pdxdt+ I |ug - ?/>(-, 0) dxdt 



+ / V/ {F,{u',K')-ed^\^^{u')-^,{K')\)d,^dxdt>0. (72) 

It follows from Proposition 2.3 that 

ed, \^^{u') - <f,{K')\ ^ in Ll^{n, x (0,T)) as e ^ 0. 

Letting e tend to in (72) provides directly the fourth point (18) the fourth point in Definition 1.2. □ 

We state now the main result, which is in fact the extension of Proposition 3.3 to a larger class of 
initial data. 

Theorem 3.4 (main result) Let uo G L°°{M.), < wo < 1, and let {uq)^ a family of approximation 
of uq given by Lemma 2.1, and let (m^)^ be the induced sequence of bounded flux solution of (V^ ). Then 
under Assumption (12), 

\uau' in Lf^^{R x [0, T]), Vp G [1, oo) 
where u is the unique entropy solution to (10)- (11) associated to initial data uq. 

Proof: Let uq G L°°(R), < uq < 1, and let > 0. There exists uq.z. in BV{R), < uq^^ < 1 such 
that for all i? > 0, 

- uo\\Li{-R,R) < C{R)v. (73) 

If one regularizes lio.^' 

into Uq y using Lemma 2.1, and if one denotes by the associated unique 
bounded flux solution, we have seen that converges almost everywhere to Uy as e tends to 0. As 
previously, we denote by Wg the regularization of uq obtained via Lemma 2.1, and the unique 
associated bounded flux solution. 

Let i? > 0, then we have 



/ / lu^-uldxdt < lu" ~ ul\dxdt + / / 

Jo J-R Jo J-R Jo J-R 



uf, — u^\dxdt 

-R 

R 

— u\dxdt. (74) 



l-R 

The contraction principle stated in Theorem 1.3 yields 



rR+MT 

\u^{x,t) — u{x,t)\dxdt < T I \uq^^{x) — UQ{x)\dx (75) 



-R J-R-MT 



where M > max^ Lip{fi). 
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We denote by ({x, t) = min(l, (i? + 1 + M{T - t) - |a;|)+). One has C = 1 on {-R, R),C>0 and 
C e L^{M.) thus 

[ [ lu^'-ulldxdt < [ [ \u^-ul\Cdxdt. (76) 

Jo J~R Jo Jr 

It follows from Proposition 1.6 that for aU ip e ^"^'^(IR x (0,r)) with ip > a.c. and 'ip{-,T) ^ 0, one 
has 

T p 

u'^ {x,t) — u1{x,t)\dtip{x,t)dxdt + j \uq{x) ~ UQ ,^{x)\ip{x,0)dx 



^ / sign{u%x,t) -ul{x,t)){fi{u^){x,t) - fi{ul){x,t))dxip{x,t)dxdt 
-I J2 [ £d,\^riu')ix,t)-ip,{ul)ix,t)\d^^p{x,t)dxdt>0. (77) 

We denote by 

Ai{t) = [-R-l- M{T -t),-R- M{T -t)], 
Mt) = [R + M(T-t),R + l + AI{T-t)], 

^^^> \ -I if x>0, 

so that 

dxC{x,t) ■d{x)xA,{t){x), dtCix,t) = ^ XMwix)- 

1=1,2 1=1,2 

Taking i^{x,t) = (T - t)C{x,t) in (77) yields : 

- [ [ \u%x,t) - ul{x,t)\C{x,t)dxdt + T [ \u^^{x) - ul „{x)\C{x,0)dx 
Jo Jr Jr 

-A//" (T-t))^ [ \u''{x,t)-ul{x,t)\dxdt 
+ / / sign{u'{x,t)-ulix,t)){f,{u'){x,t)-f,{ul){x,t))dxdt 

Jo »=l,2-^A.(t) 

-/ iT-t)J2 [ sd{x)dx\^^{u'){x,t)~^,{ul{x,t))\dxdt>0. 

Jo i=l,2-^'^'(*) 

Since M > max^ Lip{fi), one has 

/ (r-i)^/ sign{u%x,t)~ul{x,t)){Mu'){x,t)- Mul){x,t))dxdt 

Jo i=12-^A.(t) 



(t) 

<m/" {T~t)Y^ [ \u^x,t) - ul{x,t)\dxdt 

Jo j=1.2"'A.(t) 



and thus 



u''{x,t) - ul{x,t)\C{x,t)dxdt <T I \ul{x) ~ul^^{x)\C,{x,Q)dx 

Jr 

I f £^{x)dx\v,{u'){x,t) - ^,{ul{x,t))\dxdt. 

Jo 



(78) 
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We deduce, using (75) (76) and (78) in (74), that 

\u%x,t) ~ u{x,t)\dxdt <T i \ul{x)- ul^^{x)\C{x,Q)dx 

J-R JTS. 

(T-t)}^ / e{f{x)d,\ ipi{u''){x,t) - ipi{ul{x,t))\dxdt 

j=l,2-'A,(t) 

/R+MT j-T j-R 

\uo^v{x) — uo(x)\dx + I I \ul,{x,t) — Uv(xTt)\dxdt. (79) 
-R-MT ' JO J-R 

We can now let e tend to 0. Thanks to Proposition 2.3, we can claim that 

lim/ (T-t)y / e^ix)^^\ip^{u')ix,t)-ip^iulix,t))\dxdt = 0. 



We also deduce from Proposition 3.3 that 

lim / / \u1{x,t) — u,y{x,t)\dxdt = 0. 

Since C(a;,0) is compactly supported, it follows from Lemma 2.1 and the dominated convergence 
theorem that 

lim / lu^aix) - ul^^{x)\({x,0)dx = / \uo{x) - Uf)^^{x)\C{x,0)dx. 



Thus (79) becomes: 



rT i-R. i-R+MT 

limsup / / \u'^{x,t) — u{x,t)\dxdt < T \uo^^{x) — UQ{x)\dx 

e^O Jo J-R J-R-MT 

+ [ \uo{x) - uoAxMx,0)dx. (80) 
Jm 

The inequality (80) holds for any > 0, and letting ly tend to leads to 

lim / / \u'^ {x,t) — u{x,t)\dxdt = 0, 



then tends to u in L\^^{M. x [0,r]) as e 0. Since < < 1 a.e., Holder inequality gives the 
convergence in L^^^(M x [0, T]) for all finite p. □ 
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